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14.  ABSTRACT 


The  objectives  of  this  study  are  to  improve  low  magnitude  (concentrating  on  M2.5-5)  regional  seismic  discrimination  by  performing  a  thorough  investigation  of  earthquake  source  scaling  using  diverse,  high  quality  datasets  from 
varied  tectonic  regions.  Local  to- regional  high  frequency  discrimination  requires  an  estimate  of  how  earthquakes  scale  with  size  Watter  and  Taykir  (200t )  developed  the  Magnitude  and  Distance  Amplitude  Corrections  (MDAC) 
method  to  cmpirkatly  account  for  these  effects  through  regional  calibration.  The  accuracy  of  these  corrections  has  a  direct  impact  on  our  abitity  to  identify  clandestine  expkisions  in  the  broad  regional  areas  characterized  by  k»w 
.seismicity.  Unfortunately  our  knowledge  at  small  magnitudes  (i.e..  mb  <  -4.0)  is  poorly  resolved,  and  source  scaling  remains  a  subject  of  on-going  debate  in  the  earthquake  seismokigy  community  Recentty  there  have  been  a 
number  of  empiricat  studies  suggesting  scaling  of  micro-earthquakes  is  non  self  similar,  yet  there  are  an  equal  number  of  compelling  studies  that  would  suggest  otherwise.  It  is  not  clear  whether  different  studies  obtain  different 
results  because  they  anatyse  different  earthquakes,  or  becau.se  they  use  different  methods.  Wc  investigate  earthquake  sources  and  seating  from  different  tectonic  settings,  comparing  direct  and  coda  wave  anatysis  methods  that  both 
make  use  of  Empirical  Green's  Function  (EGF)  earthquakes  to  remove  path  effects.  We  devetap  improved  versions  of  the  two  methods,  and  then  apply  them  to  a  range  of  data  sets. 

Our  direct  wave  anatysis  uses  kicat  recordings  of  ckisety-  kreated.  wetl-correlaled  earthquakes,  and  a  multi  taper  method  to  obtain  time-domain  source-time- fund  ions  by  frequency  division.  Wfe  only  accept  an  earthquake  and  EGF 
pair  if  they  are  able  to  produce  a  clear,  time  domain  source  pulse.  Our  improved,  and  simplified  coda  wave  analysis  method  uses  spectrat  ratios  (as  for  the  direct  waves)  but  reties  on  the  averaging  nature  of  the  coda  waves  to  use  EGF 
events  that  do  not  meet  the  strict  criteria  of  similarity  required  for  the  direct  wave  analysis.  Both  coda  and  directwave  methods  are  limited  by  the  frequency  bandwidth  with  good  signal-lo  noise,  and  recording  limitations  such  as 
record  length,  clipping  and  sample  rate.  For  both  direct  and  coda  waves,  we  fit  the  spectral  ratios  and  perform  a  gridsearch  about  the  preferred  parameters  to  ensure  the  fits  are  well  constrained.  We  then  model  the  spectral  (amplitude) 
ratio  to  determine  source  parameters  from  both  direct  P  and  S  waves.  We  investigate  the  effects  of  using  different  fitting  approaches,  and  source- mode  I  variations  on  the  results. 

W'e  analyze  three  clusters  of  aftershocks  from  the  we  tt- recorded  sequence  following  the  M5  Au  Sable  Forks,  NY,  earthquake  to  obtain  some  of  the  first  accurate  source  parameters  for  small  earthquakes  in  eastern  North  America. 
Each  ctuster  contains  a  M~  2.  and  two  amt  am  M  3.  as  welt  as  smaller  aftershocks.  Using  direct  waves,  we  find  that  the  comer  frequencies  and  stress  drops  are  high  (averaging  t()0  MPa)  confirming  previous  work  suggesting  that 
intraplate  continental  earthquakes  have  higher  stress  drops  than  events  at  plate  boundaries.  W'e  also  demonstrate  that  a  scaling  breakdown  suggested  by  earlier  work  is  simply  an  artifact  of  their  more  band  limited  data.  W'e  atso  find 
that  the  ratio  of  Energy  to  Moment  is  high,  around  10-4.  W'e  estimate  source  parameters  for  the  M5  mainshock  using  similar  methods.  The  stress  drop  and  energy/moment  ratk)  for  the  mainshock  are  slightly  higher  than  for  the 
aftershocks.  W'e  then  appty  the  coda  method  to  the  same  sequence  to  investigate  consistency,  and  uncertainties. 

The  large  (M>4)  aftershocks  of  the  1994  Northridge,  CA.  earthquake  are  our  second  choice  of  sequence  for  comparison  of  the  methtids.  These  earthquakes  already  have  source  parameters  determined  using  time  domain  methods 
(Mori  et  al..  2003)  and  so  by  applying  our  spectral  coda  and  direct  wave  methods,  we  can  compare  most  approaches  in  common  use. 

We  atso  appty  one  or  both  methods  to  a  number  of  other  earthquake  sequences  inctuding  those  ot  recent  moderate-sized  earthquakes  in  Itaty,  and  the  W'ells,  Nevada,  earthquake. 
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ABSTRACT 


The  objectives  of  this  study  are  to  improve  low-magnitude  (eoneentrating  on  M2. 5-5)  regional  scismie 
discrimination  by  performing  a  thorough  investigation  of  earthquake  source  scaling  using  diverse,  high-quality 
datasets  from  varied  tectonic  regions.  Local-to-regional  high-frequency  discrimination  requires  an  estimate  of  how 
earthquakes  scale  with  size.  Walter  and  Taylor  (2001)  developed  the  Magnitude  and  Distance  Amplitude 
Corrections  (MDAC)  method  to  empirically  account  for  these  efTeets  through  regional  calibration.  The  aecuraey  of 
these  corrections  has  a  direct  impact  on  our  ability  to  identify  clandestine  explosions  in  the  broad  regional  areas 
characterized  by  low  seismicity.  Unfortunately  our  knowledge  at  small  magnitudes  (i.e.,  mb  <  — 4.0)  is  poorly 
resolved,  and  source  sealing  remains  a  subject  of  on-going  debate  in  the  earthquake  seismology  community. 

Recently  there  have  been  a  number  of  empirical  studies  suggesting  sealing  of  miero-earthquakes  is  non-self-similar, 
yet  there  are  an  equal  number  of  compelling  studies  that  would  suggest  otherw  ise  It  is  not  elear  whether  different 
studies  obtain  different  results  because  they  analyse  different  earthquakes,  or  because  they  use  different  methods. 

We  investigate  earthquake  sources  and  scaling  from  different  tectonic  settings,  comparing  direct  and  coda  wave 
analysis  methods  that  both  make  use  of  Empirical  Green's  Function  (EGF)  earthquakes  to  remove  path  effeets.  We 
develop  improved  versions  of  the  two  methods,  and  then  apply  them  to  a  range  of  data  sets. 

Our  direct  wave  analysis  uses  local  recordings  of  elosely-loeated,  well-eorrelated  earthquakes,  and  a  multi-taper 
method  to  obtain  time-domain  source-timc-functions  by  frequency  division.  We  only  accept  an  earthquake  and  EGF 
pair  if  they  are  able  to  produee  a  elear,  time-domain  source  pulse.  Our  improved,  and  simplified  coda  wave  analysis 
method  uses  spectral  ratios  (as  for  the  direct  waves)  but  relies  on  the  averaging  nature  of  the  coda  waves  to  use  EGF 
events  that  do  not  meet  the  striet  criteria  of  similarity  required  for  the  direet  wave  analysis.  Both  eoda  and  direct- 
wave  methods  are  limited  by  the  frequency  bandwidth  with  good  signal-to-noise,  and  recording  limitations  such  as 
record  length,  clipping  and  sample  rate.  For  both  direet  and  coda  waves,  wc  fit  the  speetral  ratios  and  perform  a  grid- 
search  about  the  preferred  parameters  to  ensure  the  fits  are  well  constrained.  We  then  model  the  spectral  (amplitude) 
ratio  to  determine  souree  parameters  from  both  direet  P  and  S  waves.  We  investigate  the  effects  of  using  different 
fitting  approaches,  and  souree-model  variations  on  the  results. 

We  analyze  three  clusters  of  aftershocks  from  the  well-recorded  sequence  following  the  M5  Au  Sable  Forks,  NY, 
earthquake  to  obtain  some  of  the  first  accurate  souree  parameters  for  small  earthquakes  in  eastern  North  Ameriea. 
Each  cluster  contains  a  M-  2,  and  two  contain  M~~3,  as  well  as  smaller  aftershocks.  Using  direct  waves,  wc  find  that 
the  corner  frequencies  and  stress  drops  arc  high  (averaging  100  MPa)  confirming  previous  work  suggesting  that 
intraplate  continental  earthquakes  have  higher  stress  drops  than  events  at  plate  boundaries.  We  also  demonstrate  that 
a  sealing  breakdown  suggested  by  earlier  work  is  simply  an  artifact  of  their  more  band-limited  data.  We  also  find 
that  the  ratio  of  Energy  to  Moment  is  high,  around  10~\  We  estimate  souree  parameters  for  the  M5  mainshoek  using 
similar  methods.  The  stress  drop  and  energy/moment  ratio  for  the  mainshoek  are  slightly  higher  than  for  the 
aftershocks.  We  then  apply  the  eoda  method  to  the  same  sequence  to  investigate  eonsisteney,  and  uncertainties. 

The  large  (M>4)  aftershocks  of  the  1994  Northridge,  CA,  earthquake  are  our  seeond  ehoiee  of  sequence  for 
comparison  of  the  methods.  These  earthquakes  already  have  source  parameters  determined  using  time-domain 
methods  (Mori  et  al.,  2003)  and  so  by  applying  our  speetral  eoda  and  direct  wave  methods,  we  can  compare  most 
approaches  in  common  use. 

We  also  apply  one  or  both  methods  to  a  number  of  other  earthquake  sequences  including  those  of  recent 
moderate-sized  earthquakes  in  Italy,  and  the  Wells,  Nevada,  earthquake. 
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OBJECTIVES 

The  objectives  of  this  study  are  to  improve  low-magnitude  regional  seismic  discrimination  by  performing  a  thorough 
investigation  of  earthquake  source  scaling  using  diverse,  high-quality  datasets  from  varied  tectonic  regions. 
Local-to-regional  high-frequency  discrimination  requires  an  estimate  of  how  earthquakes  scale  with  size.  Walter  and 
Taylor  (2001)  developed  the  MDAC  method  to  empirically  account  for  these  effects  through  regional  calibration. 

The  accuracy  of  these  corrections  has  a  direct  impact  on  our  ability  to  identify  clandestine  explosions  in  the  broad 
regional  areas  characterized  by  low  seismicity.  Unfortunately  our  knowledge  at  small  magnitudes  (i.e.,  m*,  <  -4.0)  is 
poorly  resolved,  and  source  scaling  remains  a  subject  of  on-going  debate  in  the  earthquake  seismology  community. 
Recently  there  have  been  a  number  of  empirical  studies  suggesting  scaling  of  micro-earthquakes  is  non-self-similar 
(e.g.,  Kanamori  et  al.,  1993;  Abercrombie,  1995;  Mayeda  and  Walter,  1996;  Mori  et  al.,  2003;  Stork  and  Ito,  2004; 
and  Izutani  and  Kanamori,  2001),  yet  there  are  an  equal  number  of  compelling  studies  that  would  suggest  otherwise 
(e.g.,  MeGarr,  1999;  Ide  and  Beroza,  2001;  Imanishi  et  al.,  2004;  and  Prieto  et  al.,  2004).  It  is  not  clear  whether 
different  studies  reach  different  conclusions  because  they  use  different  datasets  and  scaling  varies  with  location,  or 
because  they  use  different  methods.  Sonley  and  Abercrombie  (2006)  show  that  small  variations  in  the  commonly 
used  methods  can  lead  to  significant  differences  in  results.  Even  in  regions  that  are  well  studied,  such  as  test  sites  or 
areas  of  high  seismicity,  we  still  rely  on  empirical  scaling  relations  derived  from  studies  taken  from  half-way  around 
the  world  at  inter-plate  regions. 

In  summary,  we  address  the  following  problems: 

1 .  Do  different  studies  obtain  different  results  because  they  use  different  methods,  or  because  they  analyse 
different  data  sets?  We  will  investigate  whether  coda  and  direct  wave  methods  applied  to  the  same  datasets 
provide  the  same  scaling. 

2.  Is  scaling  dependent  upon  the  tectonic  setting?  We  will  investigate  earthquakes  from  different  tectonic 
settings  and  depth  ranges,  using  the  same  coda  and  direct  wave  methods. 

3.  There  have  been  few  studies  in  intra-plate  areas  where  seismicity  is  low  and/or  in  regions  where  a 
clandestine  test  might  occur.  The  MDAC  method  currently  assumes  earthquake  source  scaling  that  was 
derived  exclusively  from  the  western  United  States.  Can  we  extrapolate  or  transport  results  from  one  region 
to  others,  or  must  we  calibrate  to  each  specific  region?  We  will  analyse  earthquakes  from  both  interplate 
(e.g.  California)  and  intraplate  (e.g.  Eastern  North  America)  regions  to  specifically  address  this  question. 

RESEARCH  ACCOMPLISHED 

Our  approach  to  obtaining  improved  source  parameters  for  small  earthquakes  focuses  on  the  direet  and  coda  wave 
methods:  to  improve  and  investigate  them  both,  and  then  to  apply  them  to  diverse  data  sets. 

Locally  recorded,  direct  waves  from  events  have  been  least  affected  by  travel  through  the  earth,  and  so  are  thought 
to  be  the  best  candidate  for  obtaining  accurate  source  parameters.  But  there  are  only  a  limited  number  of  earthquakes 
that  are  recorded  locally,  by  sufficient  stations  to  give  good  azimuthal  coverage.  Even  fewer  of  these  have  an 
equivalently  well  recorded,  very  closely  located  smaller  earthquake  that  meets  the  stringent  criteria  required  to  be  a 
good  EGF  to  remove  path  effects.  This  EGF  method  is  the  preferred  one  for  isolating  the  source,  but  concern  about 
the  quality  of  the  EGFs  is  a  major  source  of  uncertainty  in  studies  that  use  these  methods.  In  contrast,  coda  waves 
average  radiation  from  all  directions  so  single-station  records  should  be  adequate  and  previous  work  suggests  that 
the  requirements  for  the  EGF  event  are  much  less  stringent.  It  is  thus  ideal  in  regions  with  sparse  stations  and  events 
so  that  most  events  are  only  well-recorded  by  a  single  station.  Our  approach  is  to: 

1.  Develop  an  easy  to  apply  coda  wave  spectral  ratio  method  to  obtain  source  parameters  for  large  groups  of 

earthquakes 

2.  Identify  the  mainshock  -  EGF  earthquake  pairs  that  meet  stringent  criteria  for  selecting  the  EGF,  and  obtain 

source  parameters  from  the  direct  waves  for  this  subset  of  ev  ents. 

3.  Use  the  direct  wave  results  to  confirm,  and  if  necessary  correct,  the  coda  wave  results 

4.  Apply  these  methods  to  data  sets  from  a  range  of  tectonic  environments. 

5.  Determine  the  implications  of  the  determined  source  scaling  results  for  both  coda  calibration  and  regional 

discrimination  using  MDAC  and  other  similar  means  of  source  and  path-corrected  discriminants. 

So  far  we  have  developed  improved  coda  wave  and  direct  wave  methods  that  we  will  apply  consistently  to  all 
the  data  sets.  The  new  spectral-ratio  method  developed  in  this  study  makes  the  coda  method  much  easier  and  simpler 
to  apply,  greatly  extending  the  number  of  earthquakes  that  can  be  studied.  The  coda  wave  approach  was  published 
by  Mayeda  et  al,  (2007).  We  have  applied  the  coda  wave  methods  to  several  aftershock  sequences  in  Italy 
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(Malagnini  and  Mayeda  2008,  Malganini  et  al.,  2008),  and  have  begun  to  apply  it  in  eastern  North  America.  In  the 
direct  wave  study,  we  have  investigated  the  effects  of  the  various  analysis  choices  commonly  made,  and  also 
developed  criteria  for  assessing  the  quality  of  a  particular  EGF  event.  We  have  performed  a  very  detailed  direct 
wave  EGF  analysis  of  aftershocks  of  the  Au  Sable  Forks  earthquake  (eastern  North  America),  and  obtained  some  of 
the  best  constrained  source  parameters  for  earthquakes  in  this  tectonic  setting.  The  direct  wave  analysis  method  and 
results  are  described  by  Viegas  et  al.  (2009).  In  addition  to  applying  both  methods  to  these  relatively  small  data  sets, 
we  are  all  performing  a  test  by  analyzing  the  large  (M4-5)  aftershocks  of  the  1994  Northridge,  CA,  earthquake.  Mori 
et  al.  (2003)  performed  a  detailed  study  calculating  earthquake  source  parameters  for  these  earthquakes.  They  used 
time  domain  methods  and  EGF  analysis  to  calculate  stress  drop  and  radiated  seismic  energy  for  47  earthquakes.  We 
are  now  applying  our  two  frequency  domain  methods  to  the  same  earthquakes,  thus  enabling  us  to  investigate  the 
differences  between  three  different  methods.  The  principal  aims  of  this  comparison  work  are  to  search  for  any 
systematic  differences,  and  also  determine  which  analysis  choices  lead  to  the  largest  uncertainties. 

In  addition  to  the  possibility  of  differences  between  earthquake  source  spectra  obtained  from  direct  and  coda  waves, 
in  our  research  we  have  also  uncovered  a  number  of  other  possible  sources  of  bias  in  the  methods  traditionally  used 
by  the  eoda-wave  and  direct-wave  communities  in  their  modeling  of  the  spectra.  For  example,  most  direct  wave 
studies  Fit  the  spectra  individually,  whereas  the  eoda-wave  studies  typically  invert  a  group  of  events  fixing  the 
comer  frequency  of  each  earthquake  to  be  the  same  in  all  ratios  Also,  direct  wave  studies  typically  focus  on 
obtaining  the  parameters  for  the  larger  event  in  an  EGF  pair,  whereas  most  coda  studies  ratio  all  smaller  earthquakes 
to  one  mainshoek  in  order  to  obtain  parameters  for  the  smaller  events  in  the  pairs.  We  will  include  analysis  of  these 
different  approaches  in  our  comparisons.  Although  they  sound  small,  they  can  be  significant.  For  example, 
observations  of  non-similar  sealing  in  coda  wave  studies  often  involve  a  difference  between  one  or  two  larger 
earthquakes  and  a  group  of  smaller  ones  (c.g.,  Malagnini  and  Mayeda  2008,  Malganini  et  al.,  2008).  It  is  possible 
that  this  is  an  artifact  from  the  large  earthquakes  always  being  the  larger  ones  in  the  pairs,  and  the  others  the  smaller. 


The  Au  Sable  Forks  Earthquake  Sequence,  Eastern  North  America:  Direct  Waves 

Analysis  of  earthquakes  in  stable,  mtraplate,  low'  seismicity  regions  is  important  to  characterize  these  regions,  but  it 
is  hard  because  of  the  spareity  of  earthquakes  and  stations,  and  hence  useful  data.  Relatively  little  is  known  about 
earthquake  sources  in  intraplate  regions,  and  often  relationships  based  on  minimal  data,  or  simply  extrapolated  from 
interplate  regions  are  assumed  due  to  lack  of  local  information.  For  example,  Somerville  et  al.  (2002)  use  the 
recordings  of  only  3  moderate-sized  earthquakes  to  propose  source  sealing  relationships  for  earthquakes  in  the 
Northeastern  USA  and  other  stable  continental  regions.  Their  results  imply  that  earthquakes  in  the  Northeastern 
USA  have  relatively  high  stress  drops.  This  result  was  confirmed  by  Shi  et  al.  (1998)  who  analyzed  almost  50  small 
earthquakes  from  the  region.  Their  study  included  8  empirical  Green’s  function  pairs  (Shi  et  al.,  1996),  but  mostly 
individual  earthquakes.  They  also  found  a  decrease  in  stress  drop  with  seismie  moment,  at  small  magnitudes  (< 
M~3).  This  may  indicate  that  earthquake  source  sealing  is  different  in  this  intraplate  region,  or  else  may  represent  a 
limit  to  the  resolution  as  was  found  earlier  in  the  San  Andreas  Fault  plate  boundary  region  (e.g.,  Abercrombie  and 
Leary,  1993).  Shi  et  al.  (1998)  were  mostly  limited  by  data  availability  to  using  regional  recordings  at  single 
stations,  and  so  to  frequencies  less  than  about  25  Hz. 

On  April  20,  2002  an  earthquake  of  magnitude  ML  5.3  occurred  in  the  northeastern  Adirondack  Mountains  -  the  first 
earthquake  of  such  size  to  be  recorded  by  modern  regional  broadband  networks.  Following  the  mainshoek,  the 
Lamont-Doherty  Earth  Observatory  deployed  a  small  portable  network  to  monitor  aftershocks  (Seeber  et  al..  2002). 
Between  April  22  and  November,  2002,  74  small  aftershocks  detected  and  located  in  the  epieentral  area  of  the 
mainshoek.  These  data  (200  samples/s)  represent  the  best  recorded  earthquakes  in  the  region  to  date,  and  so  provide 
an  unprecedented  opportunity  to  investigate  source  parameters  in  this  intraplate  setting. 

We  calculate  the  source  parameters  using  two  methods  commonly  applied  to  direct  waves  recorded  at  local  stations: 
spectral  modeling  of  the  individual  three-component  P  and  S  waves  (e.g.,  Abercrombie,  1995),  and  the  Empirical 
Green’s  Function  (EGF)  method  (e.g.,  Mori  and  Frankel,  1990,  Abercrombie  and  Rice,  2005).  Both  methods  require 
data  with  a  high-frequency  content  to  work  well.  Many  studies  consider  the  EGF  method  superior  as  it  corrects  for 
all  path  and  site  effects  by  using  a  smaller,  collocated  earthquake  as  an  EGF.  The  EGF  method  cannot  be  used  to 
calculate  the  seismie  moment  of  the  earthquakes,  but  this  is  the  most  reliable  information  that  can  be  obtained  from 
the  individual  spectral  analysis.  Using  cross-correlation  we  identify  three  clusters  of  earthquakes  with  very  similar 
waveforms  that  include  the  largest  aftershocks  recorded  by  the  portable  network.  All  clusters  include  an  M>2,  and  a 
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number  of  Ml -2  earthquakes,  and  two  clusters  also  each  include  an  M^3  earthquake.  We  model  the  spectra  obtained 
by  dividing  the  spectrum  of  the  large  earthquake  by  the  smaller  one.  We  also  use  regional  recordings  (100 
samples/s)  to  analyze  the  M5  mainshock  using  the  largest  aftershock  (M3. 7)  as  the  EGF.  These  earthquakes  are 
closely  located,  and  the  deconvolution  produces  a  source  time  function,  but  clearly  have  a  different  focal 
mechanism,  and  so  do  not  meet  the  criteria  for  selecting  an  EGF  that  we  describe  below.  Unfortunately,  no 
aftershocks  well  recorded  regionally  meet  these  criteria  so  we  use  the  largest  aftershock,  and  interpret  our  results 
with  caution. 

It  is  not  clear  in  many  studies  how  close  the  EGF  events  are  to  a  perfect  Green’s  function.  They  are  typically  too 
small  to  obtain  focal  mechanisms,  and  the  location  uncertainties  are  larger  than  the  preferred  separation  between 
events.  We  select  only  events  that  are  located  within  the  uncertainties  of  the  large  event,  and  have  a  high  degree  of 
waveform  similarity,  determined  by  cross-correlation.  As  a  further  test,  we  transform  all  the  spectral  ratios  back  to 
the  time  domain,  and  we  only  use  pairs  where  we  are  able  to  resolve  a  clear  source  pulse.  If  a  source  time  function  is 
observ  able  in  this  way,  then  it  demonstrates  that  the  phase  components  of  the  spectra  are  also  very  similar.  This  test 
is  similar  to  the  investigation  performed  by  Mori  and  Frankel  (1990)  to  determine  how  closely  located  earthquakes 
pairs  must  be  for  the  EGF  method  to  work  well.  Multitaper  methods  have  long  been  preferred  to  calculate  source 
spectra  and  spectral  ratios  as  they  better  represent  the  frequency  content  of  the  waveforms  than  do  cosine  and  other 
tapers  (e.g.,  Park  et  al.,  1987).  Unfortunately,  until  recently  multitaper  codes  only  worked  with  amplitude  spectra 
and  so  could  not  be  used  to  perform  the  complex  deconvolution.  We  use  the  multitaper  approach  recently  developed 
by  Prieto  et  al.  (2008)  to  perform  the  full  complex  deconvolution  so  that  we  can  use  the  same  frequency  analysis  to 
obtain  both  the  best  spectra,  and  retrieve  the  source  time  functions. 

We  first  calculate  source  parameters  using  the  individual  spectra,  to  determine  the  seismic  moments,  and  also  to 
compare  the  results  of  this  standard  method  w  ith  that  of  the  preferred  EGF  method  Before  modeling  the  individual 
amplitude  spectra,  we  re-sample  them  on  a  logarithmic  scale,  so  that  the  fits  were  not  biased  to  the  high  frequencies. 
We  fit  the  displacement  amplitude  spectra  with  the  omega-square  source  model  with  the  sharper  comer  preferred  by 
Boatwright  (1980)  and  Abercrombie  (1995)  to  obtain  the  seismic  moment,  corner  frequency  and  attenuation  Quality 
factor  (0.  We  perform  a  grid  search  around  the  preferred  parameters  to  investigate  the  uncertainties  within  a  range 
of  variance  of  fit  of  +5%.  We  take  the  mean  of  all  available  components  and  stations  to  calculate  the  final  value. 
Following  Abercrombie  (1995)  we  assume  the  circular  fault  model  of  Madariaga  (1976)  to  calculate  source  radius 
(r),  and  the  solution  of  Eshelby  (1957)  for  a  circular  crack  to  calculate  stress  drop  (Aa).  We  calculate  the  radiated 
seismic  energy  ( E$ )  by  integrating  the  velocity-squared  spectrum,  using  the  data  within  the  available  band  with  and 
the  best  fitting  model  to  extend  the  frequency  range  (Abercrombie  1995). 

We  then  apply  the  EGF  method  to  the  earthquakes  in  the  clusters  we  identified.  We  use  the  same  amplitude  spectra 
to  calculate  spectral  ratios  between  the  different  events,  and  we  also  use  the  extension  of  the  multitaper  method  to 
the  complex  spectra  to  deconvolve  the  spectral  ratios  back  to  the  time  domain  and  so  obtain  estimates  of  the  source 
time  functions.  We  only  continued  the  analysis  if  we  obtained  a  clear  source  pulse,  justifying  our  choice  of  EGF 
event.  We  resample  the  spectral  ratios  on  a  logarithmic  scale  (in  the  same  way  as  the  individual  spectra)  and  we 
model  the  spectra  ratios  using  the  same  source  model  as  the  individual  spectra.  We  use  reasonable  constraints  on  the 
fitting  parameters  as  a  further  constraint  on  the  use  of  EGF  pairs  and  ratios.  We  then  calculate  source  parameters 
from  the  spectral  fitting  as  for  the  individual  spectra. 

The  smallest  earthquakes  have  comer  frequencies  that  are  clearly  outside  the  available  bandwidth  (>80  Hz).  The 
EGF  results  are  preferred  because  of  their  better  correction  for  path  and  site  effects.  The  individual  fitting  tends  to 
underestimate  the  comer  frequencies  (and  hence  the  stress  drops)  as  they  get  close  to  the  high  frequency  limit 
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(Figure  1).  The  high  stress  drops  we  obtain  arc  consistent  with  the  previous  studies  which  Find  relatively  high  stress 
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Figure  1:  Stress  Drop  measurements  in 
Eastern  North  America.  Comparison  of  the 
Au  Sable  sequence,  with  individual  Fitting 
and  previous  studies.  The  two  studies  shown 
in  grey  are  limited  to  30  Hz.  The  individual 
measurements  underestimate  stress  drop 
near  the  bandwidth  limit  compared  to  the 
EGF.  Arrows  indicate  minimum  estimates. 
No  evidence  for  a  breakdown  in  sealing  is 
seen. 


Seismic  Moment  (N  m) 


drops  for  intraplate  regions  such  as  eastern  North  America.  They  are  consistent  with  the  results  of  Shi  et  al.  (1998) 
for  their  EGF  events,  but  we  do  not  see  a  breakdown  in  constant  stress  drop  scaling  at  lower  magnitudes,  suggesting 
that  this  was  an  artifact  of  the  limited  bandw  idth  available  to  Shi  et  al. 


Wc  obtain  relatively  high  stress  drops  and  apparent  stress  (proportional  to  the  ratio  of  Es/Mn ),  compared  to  previous 
studies,  mostly  in  interplatc  environments  (Figure  2).  This  is  consistent  with  the  hypothesis  that  faults  in  intraplatc 
settings  with  lower  strain  rates,  lower  cumulative  deformation,  and  longer  healing  time  between  ruptures  will  be 
stronger  than  faster  moving  faults  at  plate  boundaries.  Our  results  also  show  some  gradual  decrease  in  stress  drop 
and  apparent  stress  with  decreasing  moment.  Unfortunately,  this  is  largely  dependent  on  the  uncertainties  in  the 
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Figure  2a.  Comparison  of 
the  Au  Sable  Forks 
results  with  previous 
work  mainly  from  the 
western  USA.  Figure 
updated  from 
Abercrombie  and  Leary 
(1993),  corrected  for 
different  source  models. 
Note  that  recent  EGF 
studies  (dark  blue)  all 
have  higher  stress  drop 
than  earlier  studies. 


423 


2009  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


•  This  Study 

♦  Matsuzawa  (2004) 

*  Ide  et  al  (2004) 

■  Oye  et  al  (2005) 

•  Yamada  et  al  (2007) 


★  Imamshi  &  Ellsworth  (2005) 

*  Abercrombie  &  Rice  (2005) 

▼  Venkataraman  et  al  (2006) 

▼  Venkataraman  and  Kanamori  (2004) 

♦  Mon,  Abercrombie  &  Kanarrori  (2003) 


+  Mayeda  and  VNfa Iter  (2005)  *  Ide,  Beroza,  Prejean,  &  Ellsworth  (200c 

Myy 

-2  -1  0  1  2  34  5  6  7  0 

f-3 
r  “4 
r  -5 
r  -6 

r7 

t  _ ^ _ f _ (  _ i _ ^ _ ^  _ j  _ i  _ ( _ ^ _ | _ (  _ ^ _ ^ _ (  _ L  _g 

6  "7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22 

tog  M*  (Nm) 


Figure  2b.  Comparison  of 
the  Au  Sable  Forks 
results  with  previous 
work  mainly  from  the 
western  USA.  Figure 
updated  from  Ide  and 
Beroza  (2001). 


values  for  the  two  largest  earthquakes.  These  are  less  reliable  because  of  the  lack  of  an  acceptable  EGF  event. 


The  Au  Sable  Forks  Earthquake  Sequence,  Eastern  North  America:  Coda  Waves 

After  obtaining  reliable,  high  precision  direct  wave  estimates  of  the  source  parameters  of  the  earthquakes  in  the  Au 
Sable  Forks  sequence,  we  then  applied  the  coda-wave  methods,  in  order  to  compare  them.  We  began  by  performing 
a  calibration  study  following  the  procedure  of  Mayeda  and  Walter  (1996).  For  this  we  could  only  use  regional 
recordings,  and  so  only  the  9  largest  earthquakes  that  were  recorded  at  regional  distances.  This  analysis  is  therefore 
limited  to  the  40  samples/s  and  100  samples/s  of  the  regional  stations.  To  obtain  the  calibration  spectra,  it  was 
necessary  to  use  an  apparent  stress  of  at  least  10  MPa,  significantly  higher  than  the  values  typically  used  in  previous 
studies  on  the  west  coast.  This  is  consistent  w  ith  the  direct  wave  estimate  of  relatively  high  stress  drop.  The  coda 
methods  essentially  follow  the  source  parameterization  and  model  developed  by  Walter  and  Taylor  (2001).  In  this 
approach,  the  apparent  stress  and  stress  drop  are  not  independent  parameters,  as  they  are  in  the  methods  typically 
used  for  direct  wave  analysis.  This  difference  is  more  related  to  tradition  than  to  any  requirements  of  the  different 
data.  In  the  Walter  and  Taylor  (2001),  MDAC  approach,  the  apparent  stress  is  simply  proportional  to  the  stress  drop, 
and  the  radiated  seismic  energy  is  dependent  on  the  source  model  chosen  In  direct  wave  studies,  radiated  seismic 
energy  is  calculated  relatively  independently  of  the  source  model,  and  hence  apparent  stress  is  not  directly 
proportional  to  stress  drop.  After  performing  the  calibration,  we  then  calculated  spectral  ratios  of  all  8  aftershocks  to 
the  mainshock,  following  the  methods  of  Mayeda  et  al.  (2007),  and  Malagnini  and  Mayeda  (2008),  Figure  3.  We 
compare  the  comer  frequencies  obtained  from  the  direct  wave  and  coda  wave  studies  (Figure  4).  They  agree 
relatively  well  for  the  larger  events,  but  the  coda  method  significantly  underestimates  the  comer  frequency  and 
stress  drop  for  the  smaller  events.  This  is  most  likely  a  bandwidth  limitation  since  the  direct  wave  measurements 
were  made  on  the  local  recordings  which  have  larger  signal  and  frequency  bandwidth  than  the  regional  records.  This 
result  serves  as  a  reminder  that  the  coda-wave  analyses  are  subject  to  similar  problems  and  biases  as  a  result  of 
limited  bandwidth  as  have  long  been  recognized  in  direct-wave  studies. 

The  next  step  is  clearly  to  use  the  local  records  to  calculate  coda  spectral  ratios.  However,  this  requires  investigation 
into  the  use  of  higher  frequency  records,  resulting  from  smaller  volume  scattering  in  shorter  time  windows  than  have 
previously  been  used  for  coda  analysis.  The  window  length  is  limited  by  the  triggered  recordings  of  portable 
deployments,  and  the  more  rapid  decay  of  higher  frequency  signals. 

Does  it  Matter  how  the  Spectral  Ratios  are  Modeled? 

Mayeda  and  Malagnini  (2009)  analyzed  the  source  parameters  of  the  Wells,  Nevada,  M6  earthquake  and  6  of  its 
largest  aftershocks  using  the  coda  spectral  ratio  method.  They  found  non  self-similar  scaling  when  fitting  following 
the  MDAC  procedures,  and  a  joint  inversion  as  described  by  Malagnini  and  Mayeda  (2008)  and  Malganini  et  al., 
(2008).  In  order  to  investigate  the  effects  of  using  different  fitting  approaches,  we  fit  the  same  spectra  using  the 
fitting  methods,  and  two  different  source  models  that  we  are  using  for  the  direct-wave  analyses,  the  Brunc  (1970) 
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model  used  in  MDAC,  and  the  sharper-eomered  version  of  the  same  (Boatwright  1980).  The  results  are  shown  in 
Figure  5.  The  three  different  fitting  methods  all  give  broadly  similar,  overlapping  results,  but  do  show  significant 
variation,  of  the  order  of  40%  in  seismic  moment,  50%  in  comer  frequency,  and  100%  in  stress  drop.  The  variation 
in  seismic  moment  is  perhaps  the  least  expected.  In  all  ratios,  an  amplitude  ratio  between  the  large  earthquake  (of 
known,  fixed  moment)  and  the  smaller  earthquake  is  calculated.  In  the  sharper-eomer  model,  this  long-period 
amplitude  level  is  smaller,  and  so  the  moments  of  the  aftershocks  are  larger. 


78’W  76’W  74’W  72’W  70’W 


78’W  76’W  74’W  72’W  70’W 


Figure  3.  Coda  wave  analysis  of  the  Au  Sable  Forks  earthquake  sequence,  (a)  map  of  the  earthquakes  and 
stations  used,  (b)  calibrated  source  spectra  of  the  M5  earthquake  (blue)  and  the  8  larger 
aftershocks. 
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EGF  vs  Coda  Comer  Frequency 


Figure  4.  Coda  spectral  ratios  of  the  Au  Sable  Forks  earthquake  with  the  larger  aftershocks  (left)  and 

comparison  of  the  direct  wave  and  coda  wave  estimates  of  the  corner  frequencies  for  5  earthquakes 
in  the  Au  Sable  Forks  sequence.  The  coda  data  do  not  have  sufficient  high  frequency  bandwidth  to 
resolve  the  corner  frequencies  of  the  smaller  earthquakes  well. 


Results  for  Wells  Coda  fitting 
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Figure  5.  Fitting  the  spectral  ratios  of  the 
Wells  earthquakes  using  three  different 
methods:  combined  inversion  following 
MDAC  (green),  and  individually  Fitting 
ratios  using  Brune  (1970)  model  (blue)  and 
the  sharper  corner  version  (red). 


The  1994  Northridgc,  CA,  Aftershock  Sequence  -  A  Comparison  of  Multiple  Methods. 

Mori  et  al.  (2003)  used  EGFs  for  47  large  (M4)  aftershocks  of  the  M6.7  Northridge,  CA,  earthquake  to  obtain 
well-resolved  source  time  functions.  They  then  calculated  the  stress  drop,  radiated  energy  and  apparent  stress  using 
time-domain  measurements.  They  found  constant  stress  drop  but  a  possible  increase  in  apparent  stress  with 
increasing  moment.  Mori  et  al.  systematically  searched  hundreds  of  possible  small  earthquakes  to  find  the  best  EGFs 
and  so  this  data  set  provides  an  excellent  opportunity  to  compare  multiple  methods.  We  apply  our  direct-wave 
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methods  and  coda-wave  methods  to  the  same  earthquakes  and  can  compare  them  to  one  another,  and  the 
time -domain  measurements,  without  the  need  to  hunt  for  EGFs.  As  is  usual,  the  large  aftershocks  are  the  large 
earthquakes  in  the  ratios  in  the  direct  wave  study,  whereas  in  the  coda-wave  study  these  earthquakes  are  the  smaller 
in  ratios  to  the  mainshock.  So  far  we  have  calculated  very  preliminary  results  for  a  subset  of  earthquakes,  which 
show  some  systematic  differences.  We  are  waiting  to  complete  the  study  before  we  draw  any  real  conclusions  as  to 
where  the  most  significant  uncertainties  occur. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Work  is  proceeding  on  a  number  of  fronts  simultaneously,  and  wc  look  forward  to  completing  the  various  studies 
and  drawing  the  results  together  soon.  So  far  we  have  performed  the  following  analyses: 

•  Northeastern  USA  (M 1  -5):  Direct  waves  and  preliminary  coda  waves 

•  Northndge,  CA  (M4-6):  Preliminary  direct  and  coda  waves  (existing  time  domain  results  for  comparison) 

•  Wells,  CA  (M4-6):  Coda  waves,  and  effect  of  spectral  modeling  approach 

•  Various  sequences  in  Italy  (M4-6):  Coda  waves 

In  eastern  North  America  the  direct  and  coda  wave  measurements  find  High  Act  and  High  £$,  consistent  with 
intraplate  setting.  We  also  find  no  evidence  for  a  breakdown  in  scaling.  The  Coda  and  direct  wave  measurements 
agree  well  for  large  magnitudes,  but  we  need  to  use  local  data  for  the  smaller  earthquakes.  To  date,  coda  wave 
methods  have  not  been  applied  to  high  frequency,  locally  recorded,  triggered  waveforms. 

In  Italy  and  California,  studies  of  mainshock-aftcrshock  sequences  using  coda  waves  continue  to  find  that 
earthquake  scaling  is  not  self-similar  with  decreasing  apparent  stress  with  decreasing  moment. 

The  next  steps  are  to  complete  these  ongoing  studies  and  then  compare  the  results,  to  determine  any  biases  and 
uncertainties  arising  from  the  choices  of  data  and  analysis  used.  We  will: 

•  Complete  application  of  coda  spectral  ratio  method  to  northeastern  USA  earthquakes  -  including  local 
recordings  of  smaller  earthquakes 

•  Complete  application  of  both  methods  to  Northridge  Sequence  (and  Wells?)  over  range  M  to  improve 
comparison 

•  Include  fitting  coda  spectral  ratios  with  direct  approach  and  vice  versa. 

•  Investigate  bandwidth  limitations  on  fitting  more  carefully 

•  Quantify  comparisons 

•  include  results  in  MDAC. 
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